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Abstract 

We derive a new master equation with and without Markovian approximation 
for the reduced density operator of a quantum system. We study a multiphoton 
Jaynes-Cummings model (JCM) with the dissipation at a finite temperature by 
making use of the new master equation under Markovian approximation. We present 
an exact solution of a non-resonant multiphoton JCM with the dissipstion by the 
use of the superoperator technique, and obtain analytic expressions for the atomic 
inversion and the intensity of the cavity field in the JCM. It is shown that the 
new master equation under the Markovian approximation can better describe the 
suppression of the revivals of the inversion and oscillatary behavoirs of the photon- 
number distribution in the JCM. When the damping vanishes, the usual expressions 
for the inversion and the intensity are recovered. 
PACS number(s): 42.50.Md, 42.50.Dv, 03.65.Bz 



1 Introduction 



It is well known that the Jaynes-Cummings model (JCM) [1] is a model of fundamental 
theoretical importance, as the simplest nontrival model of two coupled dissimilar quantum 
systems (the atom and the field). Without dissipation, it is exactly solvable under the 
rotating-wave approximation. The solution in the presence of dissipation is not only of 
theoretical interest, but also important from a practical point of view since the dissipation 
would be always present in any experimental realization of the model. It was found 
that the collapses and revivals of the inversion oscillations predicted by the JCM are in 
agreement with the experiments done with Rydberg atoms in a microwave cavity [2-4]. 
In these experiments, the damping of the cavity mode is not negligbly small. Thus for 
a detailed comparison with experiments, the effects of cavity losses must be taken into 
account. 

In the past few years, a number of authors have treated the JCM with the dissipation 
by the use of analytic approximations [5,6] and numerical calculations [7-11]. With the 
dissipation one has to solve a master equation for the reduced density operator of the 
system. It is well known that this master equation is more difficult to solve. To our 
knowledge, no analytic solution for this model has been given for general initial conditions 
and temperature of the environment. Agarwal and Puri [12] present an analytic solution 
for the initial state of the light field being a vaccum state. Daeubler et al. [13] found an 
analytic expression for the atomic inversion and the intensity of the field in the JCM with 
the dissipation when the reservoir is at zero temperature. 

In the present paper we derive a new master equation of the reduced density operator 
for the system with and without Markovian approximation. This master equation differs 
from the previous one [5,6,14,15] due to different coupling between the system and the 
reservoir and the presence of the renormalization term in the total Hamiltonian. It enable 
us to discuss the influence of the reservoir temperature on nonclassical effects in the JCM. 
We show that the JCM is exactly sovable for the new master equation under the Markovian 
approximation, and analytic expressions of the atomic inversion and the intensity of the 
field can be obtained for an arbitrary value of the damping within the approximations 
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used in deriving the master equation. When the damping vanishes, the usual results are 
recovered. 

This paper is organized as follows: In Sec. II we derive the new master equation for 
the reduced density operator of the system with and without Markovian approximation. 
In Sec. Ill we present an exact solution of the new master equatiuon under the Markovian 
approximation for the non-resonant multiphoton JCM where the field is initially in a 
coherent state and the atom in the excited state by the use of the superoperator technique. 
In Sec. IV we calculate the atomic inversion and the intensity of the field and discuss the 
limit of a vanishing damping. We close with a summary. 

2 Master Equation 

We consider a system described a Hamiltonian H and interacting with a heat-bath envi- 
ronment (a reservoir) which consists of an infinite set of harmonic oscillators. We assume 
that the system interacting with the environmemt can be described by the total Hamil- 
tonian 

Ht = h + y,{-&- + \ m ^) + w y. + W Y. ^ (i) 

where the second term is the Hamiltonian of the reservoir, the third one represents the 
interaction between the system and the reservoir with the coupling constant hCt, and 
the last one is the renormalization term [16]. Here we have adopted a simple coupling 
between the system and the reservoir such that it satisfies the condition [V, H] = which 
is of importance in the back-action evading and quantum-nondemolition schemes [17] and 
in the study of decoherence of quantum system [18-20]. 

For later convenience, we re write the total Hamiltoian (1) in the form 

H T = H + H R + V (2) 

where 

H R = Yl h "i b tbi, V = HY J F l (3) 

i i=l 

with 

F 1 = h'£C i b i , F 2 = hJ2Cibt, F 3 = h 2 Hj2— (4) 

i i i U i 
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where bi, bf are the boson annihilation and creation operators for the reservoir. 

The total density operator for the system and the reservoir satisfies the Liouville-von 
Neumann equation in the interaction picture [15,21] 

f = < 5 > 

It is easy to that the above equation has the formal solution 

x(t) = +^Jjyw,x(o)] m + (^) a jTdf jTd^co, [v(nx(t")}} (6) 

Tracing both sides of the above equation over the reservoir, we obtain the change in the 
reduced density operator s = Tr R \ f° r the system 

s(t)-s(to) = 4 fdt'TrRiH^m'),^)^)} 
in J to i 

+ ^b 2 f dt ' f ' M"Tr R [V(t>), [V(n X (t")}} (7) 
in Jt Jt 

where we have assumed that the system and the reservoir are uncoupled and %(0) can be 
facterized as %(0) = s(t o )f(0) at the initial time. 

We assume that the reservoir is initially in thermal equilibrium at a temperature T. 
Then it is straightforward to see that 

Wr = (F 2 ) R = 0, (F 3 ) R = Au'H (8) 

where 

a / f°° i ^HICHI 2 

Au' = ih / du K n y n 9 

Jo UJ 
where J (on) is the spectral density of the reservoir. 

The total density operator [22,23] can be expressed as 

X (t) = s(t)f(t) + Xc (t) (10) 

where Xcif) represents the correlation between the system and the reservoir described by 
the density operator f(t) at time t. The reservoir assumption allows us to take f(t) = /(0). 
Then, Eq.(7) becomes 

s(t) - s(to) = -(t - t )Au'(H[H, s(to)] + [H, s(t )]H) 

H k )2 ll dt ' jC dt "{ Tr ^ V ^ [V(ns(t")f(0)}} + Tr R [V(t>), [V(t"), xMM 
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where the second term is a higher-order term than the first one and can be neglected for 
the weak damping case [21]. After dropping the second term in Eq.(6), differentiating 
both sides of Eq.(ll) with respect to time, we have the equation 

ds(t)dt = -Au'(H[H, s(t )} + [H, s{to)]H) 

H k f iy T ^[ V (t),[V(t'),s(t')f(0)}} (12) 

where the t' integration is over the correlation functions of the reservoir, which are 
charecterized by a time which is short but finite. 

Substituting the interaction term V = J2i HF { with a change of variable t' — t — r into 
Eq.(12), under Markovian approximation s(t') = s(t), we can rewrite Eq.(12) as 

s(t) = -Au'(H[H,s(t )} + [H,s(t )]H) 

/"CO _ POO 

-[H, [H, s(t )}} £ / drWy(r) - [H, / dr[Wy(r) - W^r)] (13) 

ij 1 ■:, 

where the correlation functions of the reservoir are defined by 

Wy(r) = Tr R F i (t)F j (t-T)f(0) = {F i (t)F j (t-T)) R (14) 
Wj(r) = rr fl F i (t-r)P>(t)/(0) = (F i (t-r)F i (0>i i (15) 

where all operators are in the interaction picture. These correlation functions have been 
calculated in ref.[15]. 

It is not difficult to see that under the higher temperature approximation in which the 
temperature is assumed high enough to so that the Markovian approximation is valid, 
Eq.(13) reduces to 

s(t) = - 7 [fr, [H, s(t )}} - Au[H, s(t )]H (16) 

where 7 and Ao; are constants which depends on the temperature and the spectral density 
of the reservoir, and given by 

7 = Au/ + — lim y n v n , 17 

^W + ^fjJ*!* (18) 
where P is the Cauchy principal part of the integration [15]. 



Eq.(16) is the master equation under the Markovian approximation in the interaction 
picture for the system. It is easy to convert it to that in the Schrodinger picture with this 
form: 

Ml = i[lf , p(t)} - >y[H, [H, p(t)}} - Au[H, p{t)] (19) 

which is the desired master equation for thr reduced density operator of the system in the 
Schrodinger picture under the Markovian approximation. 

Notice that in the derivation of the master equation (14) we have used the Markovian 
approximation s(t') = s(t). Without use of the Markovian approximation, through the 
transformation of variable t' = t — r in Eq.(12), then taking the approximation s{t — r) = 
s(t) — rds(t)dt where r is short but finite. To take the result under the Markovian 
approximation (16) as ds(t)dt here, we then obtain the non-Markovian master equation 
for the reduced density operator of the system in the Schrodinger picture with the form 

^ = ^[H,p]- 1 [H,[H,p]]-Au J [H,p]H 

-ja[H, [H, [H, [H, p}}}} - Aua[H, [H, [H, p]H]] 

->yP[H, [H, [H, p]]]H - AuP[H, [H, p]]H - AcuP[H, [H, p]H]H (20) 



where 7 and Auj are given by Eq.(17) and a and (5 are defined by 

a= ^ v r A jJtwt, /> = art A [J(u) | C7((i;) |», u (2D 

ll JO UJ CLLO 



Comparison Eq.(20) with Eq.(19) it is easy to see that the second and third terms on 
the rhs of Eq.(20) are the Markovian part, and the last four terms are the non-Markovian 
contribution. 

3 Exact Solution for the JCM with Dissipation 

In this section we derive an exact solution of the master equation under the Markovian 
approximatiom for the non-resonant multiphoton JCM. If we neglect the Lamb shift term 
in Eq.(14), the master equation under the approximation becomes 

dp(t) 1 



dt ih 



[H,p(t)]-^[H,[H,p(t)}} (22) 
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It is noticed that Eq.(22) has the same form as the Milburn's equation [24] under diffusion 
approximatiom, but they come from completely different physical mechanism. The former 
originates from the dissipation while the latter is from the uncontinuous unitary evolution. 

The non-resonant multiphoton Jaynes-Cummings Hammiltonian [25] describing an 
interaction of a two-level atom with a single-mode cavity field via an m-photon process 
in the rotating-wave approximation is given by 

~ vn n 

H = huj{a + a + — <r 3 ) + -K - mu)a 3 + %\{a_a +m + a+a m ), (23) 

where uj is the frequency of the cavity field; oj is the atomic transition frequency; A 
is the atom-field coupling constant; a and a + are the field annihilation and creation 
operators, respectively; a 3 is the atomic-inversion operator and a± are the atomic "spin 
flip" operators which satisfy the relations [cr+,<7_] = 2<r 3 and [03, o±] = ±2<x±. For 
simplicity, we take h — 1 throughout the paper. 

We now start to find the exact solution for the density operator p(t) of the master 
equation (22) applied to the Hamiltonian (23). Following the approach in refs. [26,27] 
we introduce three auxiliary superoperators R, S and T which are defined through their 
actions on the density operator, respectively, 

exp( J Rr)p(t) = T, [ -Z L H k p(t)H k (24) 
k=o K - 

exp(Sr)p(t) = ex.p(-iHr)p(t)ex.p(iHr) (25) 

exp(Tr)p(t) = exp(-7rif 2 )p(t) exp(-7r J H" 2 ) (26) 

where the Hamiltonia H is given by Eq.(23). 

It is straightforward to obtain the formal solution of of the master equation (22) as 
follows: 

p{t) = exp{Rt) exp(St) exp(ft)p(0) (27) 

where p(0) is the density operator of the initial atom-field system. 

We assume that initially the field is prepared in the coherent state | z) defined by 

oo y 00 

| z) = E exp(-- | z | 2 )^= | n) = J2Qn | n) (28) 
«=o 1 vn! n=0 
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where 

Q n = exp{- l -\z\ 2 )^=- [ (29) 

and the atom was prepared in its excited state | e), so that the initial density operator of 
the atom-field system takes this form: 

P(0)=( k) ( "' °) (30) 

We divide the Hamiltonian (23) into a sum of two terms which commute with each 
other, i.e., 

H = H + Hj, [H o ,H I }=0 (31) 
In the two-dimensional atomic basis H and Hj take the following form, respectively, 

tt fn + f 
H ° = U \ n-f 

where h = at a and 5 = \{ojq — mu). 

Similarly, the square of the Hamiltonian (23) can also be expressed as a sum of two 
terms which commute with each other, 

H 2 = A + B, [A, B] = (33) 

where the representations of the operators A and B in the two-dimensional atomic basis 
take the following forms: 

X _ ( uj 2 (u + f f + X 2 a m a +m + 5 2 \ , . 

{ ' u 2 (h - f ) 2 + \ 2 a +m a m + 5 2 ) [ > 

B - ZuJ \\(h-^)a+ m S(h-f)J 
It can be checked that operators A and B also commute with H D and Hj by the use 
of the following formulae: 

n\ (n — my. 

f(n)a m = a m f{n - m) , a m f{n) = f{h + m)a m (37) 

f(n)a +m = a +m f{n + m), a +rn f(n) = f(n - m)a +m (38) 



where f(h) is an operator function on the number operator. 

For convenience, we introduce the following auxiliary "density" operator: 

p 2 (t) = exp(St) exp(f i)p(O) 

From the definition of the superoperators and the initial condition (30), we 

f>2{t) = exp(-iHit) exp(— jtB)pi(t) exp(-^tB) exp(iHit) 

where we have used 



exp(— ~/t A) 



e {-yt[uj 2 (h+f ) 2 +A 2 a m a+ m +5 2 ]} 



q e {- 7 t[^ 2 (n-f ) 2 +A 2 a+ m a m +5 2 ]} 



/ - /y i\ ( exp[-icut(n+f)} 
exp ("^ = [ exp[-i^(h-f)] / 

The operator pi(t) in Eq.(40) is defined by 

*(*)><*(*) I® I e)(e | 

where 

rn 

| #(*)) = exp{- 7 t[o; 2 (n + -) 2 + A 2 a m a +m + <5 2 ]} | ze"^) 
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If we notice that the powers of the operator B can be written as 

1(h) o \ 

A k _(h) ) 



£2k = ( A k + {h) 



( S(h + f)A k + (h) \(h+f)A k + (h)a m 
\ X(n-f)A k _(h)a +m -5(h-f)A k _(h) 

where k takes zero or an arbitrary positive integer and 

777 777 

A+(n) = (2u5)\h + -) 2 + (2\u)\h + -fa m a +m 

2 2 

A_(n) = (2u5)\h - j) 2 + (2\cu)\h - j) 2 a +m a m 
Then we can write the operator exp(— 2^tB) in the form 



where 

/ m sinhiytJ A + (n)) 

E + (n, t) = cos(ytJ^h)) - 2u5(h + -) 7/ p (50) 

i m sinhiytJ Aih)) 

E_(n,t) = cos(ytjT^h)) - 2cvS(n - -) 7 V 1 (51) 

m sinhi'itJ A + (h — m)) 

E 3+ (n, t) = -2u\{n - |) U V +K ( 52 ) 

z \ A+in — m) 



m sinhiytJ A in)) 

E 3 _(n,t) = -2uX(n- -) K * K (53) 

2 aM_ n 



Similarly, we can write the operator exp(-iHjt) in the two-dimensional atomic basis as 

e X u(-iH I t) = ( ^+(M)a m \ (54) 

CXPl lt±lt) I ZV(Ma +m I>_(n,t) ) 1 J 



where 



n ^ ^ r^rf , ^2~ m ~+ m ^l ,sin[t(5 2 + \ 2 a m a +m )} 

D + (n, t) = cos\t(5 2 + X 2 a m a +m )] - 5 L / n (55) 



D (ft , t) = ^ + Aw,! - ^-'iif.^rr 1 (^) 



V5 2 + A 2 a 



+m,Qm 



and 



„ ^ , , sm[t((5 2 + A 2 a m a +m )l . . 

D 3 +(n,t) = -t\—^====Ji (57) 



„ ^ s .,sint(5 2 + X 2 a +m a m )} . . 

D 3 _ (n, i) = -iA L / =^ 58 

For latter use, we list the following formulae 

*~ ( It!) ) = ( E-fct ) a "' ( ) ^ ( ) ^ W 
and 

- ( ttt ) = ( ttt ) a+ " a+ "' ( £<M ) = ( ) <"> 



In the derivation of Eq.(54) we have used 



. 2k _ ( (5 2 + \ 2 a m a +m ) k \ , 
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2fc+i _ / S(5 2 + \ 2 a m a +m ) k X(S 2 + \ 2 a m a +m ) k a m \ , , 

til -y A(5 2 + A 2 a+ma m )fca+m _ 6 (p + X ^a+ m a m ) k ) [b ' 

where k takes zero or an arbitrary positive integer . 
Then, from Eqs.(49) and (54) it follows that 

exp(-^)exp(- 7 tS)=^ F3 _ ( ^ t)a+m F _^ f) j (63) 

where 

F + (h, t) = E + (n, t)D+(n, t) + E 3+ (h, t)D 3+ (n, t)a rn a +m (64) 

F_(n,t) = E_(h,t)D_(n,t) + E 3 ^(h,t)D 3 _(n,t)a +m a m (65) 

F 3+ (h,t) = E 3+ (h,t)D + {n,t) + E + {h,t)D 3+ {n,t) (66) 

F 3 _(n,t) = E 3 _(n,t)D4h,t) + E_{h,t)D 3 _{h,t) (67) 

Substituting Eqs.(43) and (63) into Eq.(40), we can obtain an explicit expression for 
the operator p 2 (t) as follows: 

*® -(£!?«!) 

where we have used the following symbol: 

%(t) =1 *i(*)X*i(*) I (*,J = 1,2) (69) 

with 

| * 1 (t)) = F + (h,t) | *(*)), | tf 2 (t) = F 3 _(n,t)a +m | tf(t)) (70) 

where | *&(t)) is given by equation (44). 

Taking into account the definition of the superoperator R, we can obtain the action 
of the operator exp(Rt) on the "density" operator p 2 (t) as follows: 

OO 1 

p(t) = Y,(^t) k -H k p 2 (t)H k (71) 

k=Q Kl 

where the operator H and p 2 (t) are given by Eqs.(23) and (68), respectively. 

In fact, up to now we have found the exact solution of the master equation (22) for the 
non- resonant multiphoton Jaynes-Cummings Hamiltonian (23) in the operator form (71). 
However, in most cases of practical interest, one needs to know the matrix elements of the 
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density operator p(t) in the two-dimensional atomic basis to calculate expectation values 
of observables. Although the form of the solution (71) is pleasant, it is unconvenient in 
use. Therefore in what follows we evaluate matrix elements of the density operator in the 
two-dimensional atomic basis. 

Since H D commutes with Hj, from equation (31) we obtain 

H k = E ( \ ) H k ~ l H\ (72) 

From Eq.(33) it follows that when I takes an even number, we have 

frk-ifri _ 1 ( (n + f )(S 2 + A 2 a m a+™)^ o \ 

^"2 \ " (n-f)(5 2 + A 2 a +m a m )'/ 2 J 1 J 

and when Z is an odd number, we have 

7 ~ 2^ ^ A(n - f )(5 2 + A 2 a+ m a m )('- 1 )/ 2 a +m -5{n - f )(5 2 + A 2 a +m a m )('- 1 )/ 2 ) {7 } 
Then, 

ihe} 1 2 V A(n)+/*(n) J (75) 



E ( J = I 

l£odd 



8[a k + (n)-a k _(n)} A^^gS 



A 



where these operators a±(n) and /3±(n) are defined by, respectively, 



(76) 



a±(h) = u(n + ^) ± V5 2 + X 2 a m a +m , (3±(h) = u(n - ± V5 2 + A 2 a+ m a ri 



(77) 



Substitution of Eqs.(75) and (77) into Eq.(72) yields that 



(78) 



where 



d\n) = \[a k + (h)±a k _(n)}, ^(n) = l -[(3 k + {n) ± (3 k _(n)} 



(79) 



For convenience, we introduce 



,(*)/ 



<.(*)/ 



(80) 
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Then, Eq.(78) becomes 

( f (k) (fi) A . a m 

Making use of equations (68) and (81), we can obtain the matrix Jvi^ k \t) = H k p 2 {t)H k 
where the matrix elements are explicitly given by 

M { S(t) = f^\n)i> 11 {t)^\n) + a m ^ ) \n)i> 21 (t)f { f ) (n) 

+f { +\n)i> l2 {t)^ ) \h)a +m + a m ^ ) \n)^ 22 (t)^ ) '(n)a +m (82) 
M ( 22 \t) = (f> ( j! ) \h)a^u(t)a m ^ 

+^^h)a +m i> l2 {m,t)f^\h) + f { *\n)i> 22 {t)f { *\h) (83) 
M^it) = (M ( S(t))+ 

= <P < L kY (n)a +m ^ ll (t)fi k \n) + f W (n)^ 21 (m,t)jf(n) 

+0 W '(n)a +m ^ 12 (m,t) ( /) (fc) '(n)a +m + /i fc) (n)^ 22 (m, t)^ Y (n)a +m (84) 



with 



<j>W(n) = X --^\n) (85) 



From Eqs.(71) and (82) to (85) we finally can arrive at the explicit expression of the exact 
solution of the master equation (22) for the non-resonant multiphoton JCM Hamiltonian 
(23). With these results, one can further evaluate mean values of operators of interest. 
In the next section, we will use it to obtain analytic expressions of the atomic inversion, 
distribution of the photon number and the intensity of the field. 

4 Atomic Inversion and Intensity of the Field 

In this section, we derive analytic expressions of the atomic inversion and the intensity of 
the field, and sdudy the influence of the dissipation on nonclassical effects in the JCM. 

The atomic inversion is defined as the probability of the atom being the excited state 
minus the probability of being the ground state, that is 

(a 3 (t)) = Tr[p A (t)a 3 ] (86) 
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where pA{t) is the reduced desity operator of the atom, and can be obtained through 
tracing over the field part in (86). 

Making use of the solution (86), it is easy to rewrite the inversion (87) as 

<*»(*)> = E ^-[(n\MfM - (n\M&n)] (87) 

k,n=Q K - 

It is not an easy task to calculate the two expectation values on the rhs of Eq.(88). From 
Eqs.(82) and (83) we get that 

{n\M { ${t)\n) = {$\n)) 2 \i> 1 {n,t)\ 2 + (g (h \n + m)) 2 \ ip 2 (n + m,t) \ 2 

+2Re{ff\n)g {k \n + m)^*(n, t)ip 2 (n + m,t)} (88) 

{n\M ( ${t)\n) = (g^\n)) 2 \Mn-m,t)\ 2 +(fi\n)) 2 \Mn,t)\ 2 

+2Re{f (k] \n)g {k) (n)ijj2(n, t)^{{n -m,t)} (89) 

where we have introduced the symbols: 



^i(n,t) = {n|*i(t)) (i = 1,2) (91) 

Here \^i{t)) are given by Eq.(70). These functions f ( ±\n) and <p (kY (n) can be obtained 
through replacing the number operator n in the operators f±\n) and 0*? (n) by the 
number n. 

Through a lengthy calculation Eqs.(89)and (90) can be written explicitly as 

(n\M { S{t)\n) = j{a 2 + k (n)[(l + 6)\F + (n)\ 2 + a 2 (n) ( n+ ^ 1 \ F ^ n + m )f 



+2Re(Fl(n)F 3 _(n + m))a{n)\ 



(n + m)\. 



+2a k Jn)a k _ (n)[(l - S 2 )\F + (n)\ 2 - a 2 (n) ( " + ™^ \F 3 _{n + m)\ 2 



+2Re(Fl(n)F 3 _(n + m))5a{n) 



(n + m)\. 



+a 2k (n)[(l - S 2 )\F + (n)\ 2 + a 2 (n) ^ + \F 3 _{n + m)\ 2 



-2Re(F;(n - m)F 3 _(w))(l + S)a(n)J {n + m)! ]}|fr(n, t)| 2 (92) 
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(n\M { §(t)\n) = ]{pl\n)[a 2 (n-m)\F + (n-m)\ 2 + (l-5) \ - j F 3 _(n)| 
4 (n — my. 



+2Re(F*(n - m)F 3 _(n))a(n -m)(l-6)< 



ni 



(n — m)\' 



+2(3 k + (n)f3 k _(n)[a 2 (n - m)\F+(n - m)| 2 - (1 - 5 2 )|F 3 _(n)| 2 



ni 



+2Re(F;(n-m)F 3 _(n))5a(n) ] l {v _ ^ 
+(3 2k (n)[a 2 (n)\F + (n - m)\ 2 + (1 + 5 2 )\F^(n)\ 2 



-2Re(Fl(n - m)F 3 _(n))(l + S)a(n)J ]}[W(w - m, t)| 2 (93) 

where the functions a±, f3±,F(n) and F 3 _(n) are given by replacing the number operator 
n in the their corresponding operators by the number n, and we have introduced 

|*(M| 2 = |<n|*(t)>| 2 

= \Qn\ 2 exp{- 7 t[c 2 (n + ^) 2 + 5 2 + A 2 ^^±^]} (94) 



a(n) 



r?'A 2 

(95) 



\ (n — m)!5 2 + n!A 2 

Substituting Eqs.(93) and (94) into Eq.(88) and taking into account the relation 
a±(ri)=/3±(n + m) we find that 

1 oo 

<*s(*)> = I El^(^ i )| 2 {I^H| 2 [(l + ^ 2 -« 2 W)exp(2 7 to 2 f (n)) 

^ n=0 

+ (2 - 2S 2 + a 2 (n)) exp(2 7 ta+(n)a_(n)) + ((i - 5) 2 - a 2 (n)) exp(2 7 ta 2 (n))] 
+ |F 3 _(n + m)| 2 [(a 2 (n) - (1 - g)2) ( w + m ) ! e xp(2 7 to 2 ,(n)) 

T7-. 

-2(1 - <5 2 + ^_^a 2 (n)) exp(2 7 to + (n)a_(n)) 



+a(n)(l + 5+ )J ^ n ^ ! a 2 (n))exp(2 7 to 2 _(n))] 



+2i?e(F;(n)F 3 _(n + m))a{n)8\ 



(n + m)\ 



ni 



[exp(2 7 to;(n)) + exp(2 7 tat (n))]} (96) 

As is well known, oscillations of the photon-number distribution in the JCM is a kind 
of nonclassical effects of the cavity field. To see the influence of the dissipation on this kind 
of nonclassical behavior, in what follows we discuss photon statistics in the radiation field 
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in the JCM. The reduced density operator of the cavity field can be obtained by taking the 
trace of the total density operator p(t) over the atomic states, that is, pF = Tr^p(t). Then, 
the probability of finding n photons in the radiation field is found to be 

= E ^^[(n\M { $(t)\n) + (n\M&\t)\n)] (97) 

n=0 K - 

where the two expectation values on the rhs of Eq.(98) have been given explicitly in 
Eqs.(93) and (94). 

Making use of Eq.(98), one can calculate easily the intensity of the cavity field. The 
result is 

1 oo 

(n) = -J2n\^(n,t)\ 2 {\F4n)\ 2 [(l + 5 2 + a 2 (n))eM^t a l(n)) 

4 n=0 

+ (2 - 25 2 - a 2 (n)) exp(2 7 ta+(n)a„(n)) + ((1 - S) 2 + a 2 (n)) exp(2 7 ta 2 (n))\ 
+ |F 3 _(n + m)\ 2 [(a 2 (n) + (1 - 5) 2 ) ^ + , m) ' exp(2 7 ta 2 (n)) 

_2(1 - S 2 + (n+ , m) ' a 2 (n)) exp(2 7 ta+(n)«_(n)) 
n! 

-a(n)(l + 5 - W (n+ , m)! a 2 (n)) exp(2 7 ta 2 (n))l 
V n! 

+2^e(F;(n)F 3 „(n + m))a(n) J ^ + , m) ' [(2 - 5) exp(2 7 to 2 h (n)) 

V n! 

-45exp(2 7 ta + (n)a_(n)) - 2 exp(2 7 to 2 (n))]} (98) 

In order to better understand the influence of the dissipation on nonclaasical effects 
in the JCM, we consider the resonant case, i.e., 5 = 0. In this case, the expressions of the 
inversion , photon-number distribution and the intensity can simplified significantly. The 
inversion reduces to 

(*»(*)) = E I Qn I exp[-4 7 AV —L] cos[2Ati A —L] (99) 

n=0 U - » n - 

The photon-number distribution becomes 

/ x i9 I , r >9 (ra + ra)!, r „, /(n + m)!,| 

p(M) = -|Q»| 2 <h + exp[-4 7 A V - ; ]cos[2Aty l - ; ][ 

1 f n 1 / n! 1 

+ -\Q n -m\ 2 1 - exp[-4 7 A 2 t- — -] cos[2AtW- — -] (100) 

2 [ (n — my y [n — my J 

And the intensity of the cavity field is given by 

(n(t)}=n + ^-^e- n J2-exp[-^X 2 t y ; ]cos[2Ati/ 1 ; ] (101) 

Z Z n=0 ' 
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where n = \z\ 2 is the initial mean photon number in the field. 

From the above expressions we see that the damping term in Eq.(22) leads to the 
appearance of the decay factors exp{— A^X 2 ^ 11 ^ 1 ^ } and exp{—A r yX 2 t ^; m y } in Eqs.(lOO), 
(101) and (102), which are responsible for the destruction of revivals of the atomic in- 
version and the weakening of oscillatary behaviors of the photon-number distribution. In 
Fig.l we plot the time evolution of the inversion. With increasing the damping 7, one can 
observe rapid deterioration of the revivals of the inversion. It is worth mentioning that 
the inversion, oscillations of photon-number distribution and the intensity are weakened 
with increasing temperature of the environment due to the relation (17). 
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FIG.l: The atomic inversion (<fs) as a function of t for m = 1, \z\ 2 = 20 and 
(a) 7 = 0.0001, (b) 7 = 0.0005, (c) 7 = 0.001. Here we have let A = 1 

It is obvious that when the damping vanishes, Eqs.(100), (101) and (102) take the 
form 

<*s(*)> = E I I 2 cos[2aJ5^S (102) 

n=0 » n - 



p(n,f) = |Q„| 2 cos 2 [At^±^] + |Q„- m | 2 sin 2 [A^^±^] (103) 
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and 



(n(t))=n + - 




(n — m)\ 



(104) 



which are just the usual expressions without dissipation. 

5 Summary 

we have derived a new master equation with and without the Markovian approximatiom 
for the reduced density operator of a quantim system. Making use of the new master 
equation under the Markovian approximation, we we have studied the JCM with the 
dissipation at a finite temperature. We have found an analytic solution of the non-resonant 
multiphton JCM with the dissipation and obtained analytic expressions of the atomic 
inversion, the photon-number distribution and the intensity of the cavity field,which can 
reduce to usual expressions when the damping vanishes. We have shown that the damping 
suppresses the revivals of the atomic inversion and the oscillatary behaviors of the photon- 
number distribution. It must be mentioned that we do not discuss applications of the 
non-Markovian equation in the present paper, but it may be useful for the study of the 
non-Markovian relaxation processes in transient optical phenomena [28]. 

We would like to thank our colleagues at Theoretical Physics Division , Nankai In- 
stitute of Mathmatics for their useful discussions and comments. This research was sup- 
ported by the National Natural Science Foundation of China. 
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